Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 12 January 2013 



(MN M£X style file v2.2) 



The role of binding energies of neutron stars on the 
accretion driven evolution 



o 

IX, 



00 

43 
6 

o3 



Manjari Bagchi 1 



1 Physics Department, West Virginia University, Morgantown, WV 26506 



12 January 2013 



> 

On 

cn . 1 introduction 

o ■ 

I ■ Millisecond radio pulsars (MSRPs) are believed to de- 
• • scend from the neutron stars in low mass X-ray binaries 
m ^ . (LMXBs) that have been spun up by acquiring angular 
momentu m through accretion of matter from their com- 
panio ns (|Alpar et alj 1 19821 ; iRadhakrishnan fc Srinivasanl 
1982). kHz quasiperiodic oscillations and millisecond 
pulse periods of n eutron stars in several L MXBs sup- 
port this theory (jWiinands fc van der Klis! (|l998l ) and 
references therein) and the most compelling proof came 
from the MSRP J1023+0038 which h as lost its ac- 
cretio n disc within the past decade (|Archibald et al.l 
2009). In spite of conti nuing efforts to understand the 
evolu tion of LMXBs dWebbink. Rappaport fc Sayonik] 
Rappaport. Verbunt. fc Jossl 



ABSTRACT 

Millisecond pulsars arc believed to descend from low mass x-ray binaries. Observable 
parameters of binary millisecond pulsars e.g. mass of the pulsar, mass of the com- 
panion, spin period of the pulsar, orbital period, orbital eccentricity etc are used to 
probe the past accretion history of the millisecond pulsars. But unfortunately in these 
studies, binding energy of the neutron star is not considered commonly. We show that 
the effect of the binding energy is significant in the estimation of the amount of mass 
accretion and thus should be incorporated in models for binary evolutions. Moreover, 
different Equations of State for dense matter give different values for the accreted 
mass for the same amount of increase in the gravitational mass of the neutron star 
implying the need of constraining dense matter Equations of State even to understand 
the spin-up procedure properly. 
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iBelczvnski et all 
2008), sometimes discoveries of MSRPs with peculiar 



ILin et all [2010). But all these evolutionary models usually 
ignore the effects of binding energies of neutron stars which 
should be taken care of for a more realistic understanding 
of the binary evolution as neutron stars posses significant 
amounts of binding energies. In the present work, we show 
that the change in the total binding energy in accretion 
driven evolution is significant and should be incorporated 
in realistic models of binary evolutions. 



2 STUDY OF BINDING ENERGY OF 
NEUTRON STARS 

The gravitational mass of a neutron star is defined 
as Mg = J 47rr 2 p(r) dr where p(r) is the mass 
density at a radial distance r from the center and 
R is the radius of the neutron star. This is the 
measurable entity and being u sed in studies of neu- 
tron star mass distributions ( Thorsett & Chakr abortvl 
19991: I Schwab, Podsiadlowski, & Rappaort 
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characteristics impose problems and compel us to think 
of alternative scenarios. Two such systems are binary 
MSRPs PSR J1 9 03+q 327 which has an eccentric orbit 
|Champion et al.1 120031 ) and PSR J1614-2230 which is 
highly recycled with c omparatively larger companion mass 
l|Demorest et al.ll2010l ): and special evo lutionary channels 
have been proposed for these objects (jFreire et al.1 l20ld: 



Kiziltan. Kottas. fc Thorsett! 2010l ). But in the case 
of relativistic objects like neutron stars, one should 
think about the proper volume element dV which be- 
comes 4nr 2 [l — 2Gm G (r)/Vc 2 ] ^ 2 dr in Schwarzschild 
geometry. So the proper mass of the neutron star be- 
comes Mp = J Q fl 4-Kr 2 [l — 2Gm G (r)/rc 2 ] 1//2 p(r) dr 
and the total number of particles (baryons) be- 
comes Nb = Xf" 4-7rr 2 [l — 2Gm G (r)/rc 2 ] n(r) dr 
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where n(r) is the particle number density at 
r. This makes one define the baryonic mass as 



Af s = m B J 47rr 2 [1 - 2Gm G (r)/rc 2 



' n(r) dr where 



mp is the mass of a baryon. The binding energy of the 
neutron star is BE — (Mp — Mg) c 2 where c is the velocity 
of light and the gravitational binding energy of the neutron 
star is BEg = (Mp — Mg) c 2 . As the only difference in 
the expression of Mp from that of Mq is the appearance 
of the metric element [l — 2Gm G (r)/rc 2 ] inside the 
integral, so BEg is actually the manifestation of the space- 
time curvature for objects as compact as neutron stars. 
Solving Tollman-Oppenheimer-Volkov (TOV) equations of 
hydrostatic equilibrium, one can obtain the values of all 
these entities with a suitable choice of the Equation of 
State (EoS). It is noteworthy to mention here that although 
TOV equations have been derived for non-rotating stars in 
Schwarzschild geometry, these are good enough even for the 
millisecond pulsars. The definition of the proper volume dV 
is different for a rotating-star where the metric is different, 
but we have found with the freely available "Rapidly 
Rotating Neutron Star" code by Nikolaos Stergioula^j that 
for a fixed value of Mg , Mp usually decreases by an amount 
less than 1% if we increase the angular velocity uu from 
zero to 3000 s . So in the present work, we keep ourself 
confined within Schwarzschild geometry. One important 
question that arises here i s wha t is the value of mp ? 
Following I Arnett fc Bowers! (I1977T ). we set m B equal to one 
atomic mass unit "«" (based on C 12 scale, 931.49402 MeV) 
which is appropriate as pre-collapse evolved stars contain 
large amounts of C 12 or nuclei of similar binding energy per 
nucleon. The mass of H 1 is 1.007825 u and the mass of He 4 
is 4.002602 u, so choosing m B = 1 u is valid even for the 
present work although we are mainly concerned about the 
change in Mp due to the accretion of hydrogen or helium 
rich matter on the neutron star from its binary companion 
(usually in the red-giant phase). 

Among nu merous nuclear EsoS, we choose a stan - 
dard EoS APR (|Akmal. Pandharipande. fc Ravenhaiilll998l ) 
which agrees with the possible Mg — R regions for three neu- 
tron stars det ermined from X-ray observat ions of thermonu- 
clear bursts (|6zel. Bavm. fc Giiverl |2010| ) as well as gives 
maximum value of Mg to be greater than the mass of the 
most massive neutron star discovered till date - the millisec- 
ond pulsar PSR J1614 -2230 having M G = (1.97 ±0.04) M© 
l|Demorest et al.ll2010[) . But there is a possibility of discover- 
ing even heavier neutron stars in future - the "black widow" 
pulsar B1957+20 may be one such example which has the 
lower limit of M g as 1.66 Mq with the best fit value a s 
2.40 ± 0.12 Mq l|van Kerkwiik. Breton, fc Kulkarnil I20T0I) . 
Further study is needed to get a precise mass measure- 
ment for this neutron star. But at the present situation, 
one should not exclude very stiff EsoS in any study involv- 
ing EsoS. As an exam ple of such a stiff EoS, we use MS0 
(Miillcr & Scrot 199 jj) although it disagrees with the find- 
ings of Ozel. Bavm. fc Giiverl (|2010l ) . As an EoS much softer 
that APR would not give Mg as high as 1.97 Mq, we do 
not use any such EoS. It is worthwhile to mention he r e tha t 
the statistical method of lSteiner. Lattimer. fc Brownl (|2010T ) 
favors EsoS which are even softer than APR at high densi- 
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ties (see their figure 8). But their results fail to give Mg 
as high as 2.4 Mq and there is enough scope to improve 
their method (see their discussions) . On the other hand, the 
st udy of the cooling phase o f the X-ray burster 4U 1724-307 
bv lSuleimanov et all l|2010l ) favors an EoS as stiff as MS0. 

In the top panel of Fig.[T] we show mass-radius relations 
of neutron stars as obtained by solving TOV equations with 
these two EsoS. It is clear that the difference between three 
types of masses increases with the increase of Mg leading to 
the increase in BEg and BE with Mg (shown in the middle 
and bottom panels of the same figure). This happens be- 
cause with the increase of Mg, neutron stars become more 
compact and the general relativistic effects more dominant. 
For a fixed value of Mg, the values of BEg and BE are 
higher for a softer EoS than those for a stiffer EoS due to 
the higher compactness of the neutron star in the first case. 
This f act had been observed earlier by iPods iadlowski et al.l 
(|2005h . While BEg is useful to understand general relativis- 
tic effects, it is BE which is relevant in case of the accretion 
driven evolution which we discuss next. 

We study the role of the binding energy in the evo- 
lution of millisecond pulsars which are the results of the 
mass accretion onto the initial slow neutron stars from their 
companions. Let us assume a millisecond pulsar with grav- 
itational mass Mgj and binding energy BEf initially was 
a slow neutron star with the gravitational mass Mg,; and 
binding energy BEi before the accretion process started. 
Here Mgj > Mg,z implying BEf > BEi (bottom panel of 
Fig. [TJ. The increase in the gravitational mass is AMg = 
Mgj — Ma,i and the increase in the binding energy is 
ABE = BEf — BEi and the total amount of mass accreted is 
given by M acc = AM G + ABE/c 2 = AM p. Unfortunately, 
people usually assume M aC c = AMg; one of the main aims 
of the present paper is to point out this error. 

In Fig. [2] we plot Mg,; along X-axis, Mgj along Y- 
axis and M acc as the colour palette. The increase of Mg of 
a neutron star with the amount of accreted matter can be 
studied along a straight line parallel to the Y-axis. The X- 
coordinate of the line will give the initial value of the Mg 
(Mg,;) and the Y-coordinate at any point in the line will give 
the final value of Mg (Mgj) for which the amount of mass 
accreted can be obtained from the colour at that point. But 
Mg,; is not the observable entity, it is Mgj which can be 
measured. By noting the colours at different points along a 
straight line parallel to the X-axis, one can get the amount 
of accreted mass, where the Y-coordinate of the straight 
line should be set equal to the measured value of Mgj and 
the X-coordinate of any point along the line will correspond 
to the assumed value of Mg,;. Let us substantiate this for 
the case of the millisecond pu lsar PSR J1614-2230 which 
has M G j = (1.97 ± 0.04) M© (iDemorest et alj|2010h . spin 
period 3.15 ms and a companion of mass 0.5 Mq. Table [T] 
shows the values of M acc for different values of Mg,; setting 
M G j = 1.97 M Q . We find that M acc (AMp) is quite high 
for Mg,; << Mgj which is difficult to explain with present 
evolutionary models. There are two possibilities to explain 
the fact, first the neutron star was probably born massive 
(Mg,; > 1.8 Mq) or the accretion process is much more 
efficient. Further study is needed to be sure of possible origin 
of PSR J1614-2230. Moreover, we see from Table [Q that for 
fixed values of Mgj and Mg,;, different EsoS give different 
values for M acc . For a softer EoS, fixed values of Mgj and 



binding energy of neutron stars 3 













M P (APR]-., 


M,,(MSO) ■ 






M B (APR) . \ 


MgCMSO) 






M G (APR)~_^^ 




\\ 






A 






Figure 1. The upper panel shows Mass- Radius curves for neutron 
stars as obtained using EoS MSO and EoS APR. The values of Mb 
have been computed using mg = 1 u. The middle panel shows 
the variations of the gravitational binding energy BEq with Mq 
and the bottom panel shows the variations of binding energy BE 
with Mq for the same set of EsoS. 



Mg,i give a larger value of Macc than that obtained for a 
stiffer EoS because of the larger compactness of the neutron 
star in the first case. Thus proper understanding of binary 
evolution requires the true knowledge of dense matter EoS. 
Similarly for any other MSRP, M acc can be estimated with 
proper choices of Ma,i and the EoS. 
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Figure 2. Plot of Mq, 
MSO and EoS APR. 



Mqj and Macc as obtained using EoS 



Table 1. M acc as obtained using EsoS APR and MSO for different 
values of Mgj when Mq j = 1.97 Mq which gives Mb f = 
2.322 Mq for EoS APR and M B j = 2.237 Mq for EoS MSO. 
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3 DISCUSSIONS 

The extra contribution of the binding energy to the esti- 
mate of total accreted mass M acc makes the recycling pro- 
cess more efficient as the conservation of angular momentum 
gives If ujf = IiUi + M a cc Ra Va- Here the angular momen- 
tum of the accreted material i s given by J acc = M aC c Ra Va 
(IKonar fc Bhattacharvalll999f) . Ra is the Alfven radius and 



Va is the Keplerian velocity at Ra; Jj = If^f, Ji = h^>i 
where Jf, If, ujf are the final angular momentum, moment 
of inertia and angular velocity and Ji, U, uji are the initial 
angular momentum, moment of inertia and angular velocity. 
Moment of inertia of a neutron star increases almost linearly 
with the increase of Mg except very close to the maximum 
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mass (|Bagchill2O10l ). So if one perform the calculation con- 
sidering Mace = AM G + ABE/c 2 instead of M acc = AM G 
with the same choice of all other parameters, uij will be 
higher in the earlier case. The effect of binding energy will 
be important i n any study for P a j — Ma,f correlation as re- 
cently done bv lZhang et al. (2010) without considering this 
effect. 



4 SUMMARY AND CONCLUSIONS 

We find that neutron stars posses significant amount of 
binding energies which increase with the increase of the 
gravitational mass of the neutron stars. So when a neutron 
star gains mass by accreting matter from its companion, its 
binding energy also increases and the amount of accretion 
needed becomes greater than the difference in the gravita- 
tional masses before and after the accretion process. Indeed 
it is the difference between the baryonic masses before and 
after the accretion process. This fact makes it more difficult 
to explain the evolution of massive MSRPs like PSR J1614- 
2230. It seems that probably PSR J1614-2230 was born mas- 
sive [Mg,% > 1-8 Mq). If this is indeed the case, then it is 
quite possible of discovering more neutron stars with such 
high gravitational masses - even the neutron stars in the 
double neutron star systems may have such high gravita- 
tional masses whereas in the present sample of double neu- 
tron stars, the gravitational mass es are usually much smaller 
- wit h a mean of 1.35 Mq (|Kiziltan. Kottas. fc Thorsett! 

I2010I ). 

We also find that for the same amount of increase in 
the gravitational mass of a neutron star, different Equations 
of State predict different values for the increase in the total 
binding energy. This fact reveals the importance of knowing 
the true nature of the matter inside a neutron star for better 
understanding of the evolution of the neutron star through 
mass accretion from its companion. We use two sample EsoS 
in this work but the correctness of these are not firmly es- 
tablished. Measurement of the moment of inertia from the 
faster component of the doubl e pulsar system PSR J0737- 
3039 l|Lattimer fc Schut zH2005T ) or detection and analysis of 
low-f requency 600 Hz) gravitational waves ( Ozel et al.l 
I2010I ) might help us to constrain dense matter EsoS in near 
future. 
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